Determination of the age of the earth from Kamland measurement of geo-neutrinos 
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The low energy component of the anti-neutrino spectrum observed in the recent Kamland exper- 
iment has significant contribution from the radioactive decay of 238 U and 232 Th in the earth. By 
taking the ratio of the anti-neutrino events observed in two different energy ranges we can determine 
the present value of the Thorium by Uranium abundance ratio, independent of the U,Th distribu- 
tion in the earth. Comparing the present abundance ratio with the r-process predicted initial value 
we determine the age of the earth as a function of Am 2 and Sin 2 26. We find that the age of the 
earth determined from KamLAND data matches the age of solar system (A.bGyrs determined from 
meteorites) for the LMA-I solution. For the LM A-II solution the age of the earth does not match 
the solar system age even at 90%C.L. 

PACS numbers: 



The recent results from KamLAND p| are significant 
for establishing the LMA solution of the solar neutrino 
problem. In addition they present the first statistically 
significant measurement the anti-neutrinos from radioac- 
tive decay of 238 U and 232 Th present in the crust and 
mantle of the earth. The possibility of detecting geo 
neutrinos from radioactivity in the the earth was first 
raised by Eder and was revived by Krauss et al Q 
and Kobayashi and Fukao 0] . Quantitative spectra of the 
anti-neutrino events which could be observed at Kamioka 
and Gran Sasso from radioactivity of Th and U in the 
earth has been given by Rothschild et al Q and Ragha- 
van et al . The emphasis of these papers and the recent 
work of Fiorentini et al is to use the observations of 
v e events at KamLAND and Borexino experiments to de- 
termine the distribution of U, Th in the crust and mantle 
of the earth, and determine what fraction of the 4QTW 
heat generated by earth from radioactivity. KamLAND 
reports observing 9 v e events below 2.6MeV (visible 
positron energy) which are ascribed to U,Th radioac- 
tivity in the earth. Subsequently a analysis of the geo 
neutrinos from KamLAND observations has been done 
by Fiorentini et al with the aim of determining the 
geological distribution of U and Th in the earth. Fioren- 
tini et al. assume that [Th/U] = 3.8 the solar system 
value and try to determine the U and Th content of the 
core and mantle. 

In this paper we have a different aim which is to ex- 
tract the average thorium by uranium abundance ratio, 
[Th/U] from the KamLAND data, as this parameter is an 
ideal chronometer for measurement of cosmological time- 
scales. Th and U present in the earth are produced in 
supernova by r-process nucleosynthesis Q . The theoret- 
ical prediction initial abundance ratio is robust against 
perturbations of the astrophysical parameters at the site 
of the r-process to within ±5% 10]. By determining 
the present abundance ratio of [Th/U] in the earth from 
KamLAND, we can measure the age of the earth-i.e the 
time elapsed between the supernova explosion where the 
Th and U in the earth were produced, and the present. 



KamLAND observes geo anti-neutrinos in the (positron 
signal) energy range 0.9MeV < E < 2.5MeV. Thorium 
decay neutrinos have a maximum E value of l.bMeV 
whereas Uranium neutrinos contribute in the entire ob- 
served energy range of the geo signal. The geo neutrinos 
undergo an energy independent suppression as their av- 
erage distance from the detector ~ 10 3 fcm is much larger 
than the oscillation length ~ 10 2 km. We take the ratio 
of the geo events, in the the (0.9 — 1.5)MeV bins and the 
geo events, in the (1.5 — 2.5)MeV energy bins to deter- 
mine the [Th] /[U]. This enables us to determine the age 
of the earth independent of the geological distribution of 
Th and U in the earth. 

Using the KamLAND data we can determine the age 
of the earth as a function of Am 2 and Sin 2 29, as these 
parameters determine the reactor background that has 
to be subtracted from the KamLAND observed events 
to arrive at the geo-neutrino signal. In Fig.l we show 
the age of the earth as a function of the mixing an- 
gle for the LMA-I (Am 2 = 7 x 10" 5 e^ 2 ) and LMA-II 
(Am 2 = 15 x W- 5 eV 2 ) solutions [HI- We estimate the 
total uncertainity in the age (for a given Am 2 ,Sin 2 29) 
to be ±36.5%. The central shaded regions are the best 
fit points (x 2 ~ 0) and the two outer lines enclose the 
allowed region of parameter space with C.L. of 68.3% 
and 90% respectively. We see from Fig.l that for the 
LMA-II solution t age > 5.0Gyrs ( at 90%C.L) and 
t a ge > 9.0Gyrs (at 68.3%C.L) . It is known from the 
radiochemical dating of meteorites that the age of the 
oldest meteorites (the chondrites) is 4.5Gyrs 12]. If the 
age of the earth is close to the age of the chondrites then 
the LMA-II solution is ruled out. We see from Figl that 
for the LMA-I solution, the age of the earth agrees with 
the age of the meteorites for the mixing angle in the range 
0.88 < Sin 2 26 < 0.93. 

The KamLAND detector consists of about lkiloton of 
liquid scintillator surrounded by photo-multiplier tubes. 
Electron anti-neutrinos are detected by means of the in- 
verse beta decay v e + p — > e + + n by looking for the e + 
in delayed coincidence with the 2.2MeV 7-ray from the 
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FIG. 1: Age of the earth at 68.3C.L.% (dashed curves) and 



90%C.L. (solid curves) 
parameter space. 



Shaded region is the \ — allowed 



neutron capture by protons (n + p — > d + 7). The e + 
annihilate in the detector producing a total visible en- 
ergy E which is related to the initial v e energy,^, as 
E = E v - (m„ -m p + m e ) + 2m e = E V - 0.78MeV. Only 
those radioactive product D e 's with energies above the 
inverse beta decay reaction threshold of m n — m p + m e — 
1.8MeV can be detected. The main source of z? e 's and ra- 
diogenic heat in the earth are decays of 23S U, 232 Th and 
40 K. The v e from 40 K decay have E vmax = 1.31 MeV 
and will not register in the KamLAND detector. In the 
decay chain of 238 U, only D e 's from the (3 decays of 234 Pa 
and 21A Bi are above the threshold for detection in Kam- 
LAND. In the 232 Th decay chain, D e 's from the beta de- 
cays of 228 Ac 212 Bi contribute to the KamLAND signal. 
Thorium decay v e 's will contribute only to the (positron 
signal) energy bins below 1.5 MeV whereas Uranium v e 's 
will contribute to all energy bins below E — 2.5MeV. 
This fact enables us to separate the thorium neutrino 



signal from that of uranium. The energy spectrum of 
£/ e 's from each of these beta decays can be expressed an- 
alytically as follows, 
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where F(Z,E U ) is the Fermi function that accounts for 
the distortion of the spectrum due to Coulomb attrac- 
tion of the outgoing e~ with the nucleus, the sum is 
over each of the beta decays in the X (= Th, U) de- 
cay chain with Q value above the 1.8MeV threshold, 
A(Q, Z) are constants obtained by the normalizing the 
spectrum for each term in the sum to unity. The re- 
maining terms are kinematical factors for the two body 
decay (assuming the recoil energy of the nucleus is neg- 
ligible). Other nuclear physics effects can be parame- 
terized by adjusting the overall normalization to match 
the tabulated experimental values [13j ( we fit the nor- 
malization of U and Th spectrum by requiring that 
/ = / dEa(E)r)(E) = 0.51 (2.52) x 10~ 44 cm 2 for Tho- 
rium (Uranium) where a is the inverse beta decay cross 
section shown in © below. 

The neutrino flux from the earth at a location Rd of 
the detector can be expressed as the integral 
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where nx(f) is the number density of the radioactive 
atoms X (= U,Th) and tx is the lifetime of X and 
Pee(\r — Rd\) is the v e survival probability. Assuming 
that nx(r) is approximately constant within a spherical 
shell, the expression for the flux from a shell of constant 
density of radioactive atoms can then be written as 
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where [X], is the number of X atoms per unit mass in 
the shell i , is the mass of the i' th shell and R e is the 
radius of the earth. The geometrical factor d depends 
upon the thickness of the shell at the site of the detector 
and is given by 
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where x\ = r\/R e ,X2 = r%/R e are the inner and outer 
radii of the shell in units of R e , fi is the cosine of the angle 
between the position of the detector, Rd, and a point r 
inside the shell (the origin of the coordinates is chosen at 
the center of the earth) . The survival probability of u e is 
a function of the distance between the detector and the 
point x in the shell and is given explicitly by 



Pee = 1 - sin 2 26 
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Inserting P ee in the expression for the geometrical factor 
and carrying out the integration over the shell thickness 
x and the angular variable fi, we can express the energy 
dependence of Gi as 
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where Gi(0) depends on the shell thickness and fi(E) 
is slowly varying function of E (when Am 2 ~ 10 -4 — 
10~ 5 eV 2 ). For the continental crust, with thickness 
30fcm, G cc (0) = 3.54 and f cc {E) is a monotonically de- 
creasing function of the positron signal energy, E, with 
/ cc (0.9) = 1.00 and / cc (2.5) = 0.98. For the oceanic 
crust with thickness of about 6km, G oc = 4.34 and 
foc(E) decreases monotonically from / oc (0.9) = 0.88 to 
/ oc (2.5) = 0.82. For neutrinos coming from the mantle, 
there is an energy independent suppression, and we have 
G m (0) = 1.5 with f m (E) = 1. 

The total contribution to the neutrino flux from the 
crust and mantle can be written as 
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where the index i = cc, oc, m denotes the continental 
crust, oceanic crust and mantle respectively. In the ab- 
sence of significant chemical segregation of U and Th the 
ratio ([T/i]j/[E/]i) can be taken to be the same in the 
crust and the mantle, gi is the location parameter which 
represents the fraction of the the continental crust vis- 
a-vis the oceanic crust surrounding the detector . For 
a detector in Japan which has the oceanic crust on one 
side and the Asian continental crust on the other side we 
may take g cc ~ g oc ~ 0.5. 

The number of detection events Ni in an energy bin 
centered at i is, 
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0.9 is the threshold of the lowest energy bin). 



The ratio of the neutrino events in the energy bins be- 
tween (E min — 1.5)MeV, Nj, and the events in energy 
bins (1.5 — 2.5)M eV , N]j, depends only on the Thorium 
to Uranium ratio [T7i] /[[/], and the spectral shape of U 
and Th neutrinos folded with the cross section. The ge- 
ology factor in the square brackets in JJJ cancels out as it 
is independent of energy (to ±1%). Specifically the ratio 
of the geo-neutrino events in the two energy ranges can 
be written in a simple form 
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One can determine the ratio Nj/Njj and from that 
determine the global average abundance ratio of Th by 
(7)U using ^| From the initial r-process abundance [l(j 
[Th/U] = 1.169 ± 0.08(lfi) and the decay lifetimes of U 
and Th (tu = 6A5Gyrs, TTh = 20.03Gyrs) we can relate 
the time elapsed (in Giga-years) between between the 
supernova explosion where the earths U, Th were formed 
and the present, 



[Th/U] = [Th/U} exp[ t age ( 
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We can directly relate the age of the earth to the exper- 
imentally measured quantity Nj/Njj as 



t age = 9.45 In 
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where n p number of free target protons in the fiducial 
volume (3.46 x 10 31 for this experiment) , t is the ex- 
posure time (lA5.1days) and d e // is the detector effi- 
ciency (78.3%) and e is the width of the energy bins 
(0.425MeF). The low energy cutoff in KamLAND is 
0.9MeV. 

The cross section for the v e + p —> e + + n reaction is 
given by 01 



where E e = E v — (m n — m p ) is the positron energy and 
p e is the corresponding momentum. Thorium neutrinos 
have a E < l.hMeV, whereas Uranium neutrinos con- 
tribute in the entire range of (E min — 2.5)MeV (where 



We calculate the reactor background using the proce- 
dure given in |l5j [. We sum the neutrino flux from 16 
reactors using their power production and distances as 
input and assuming an average fuel composition in the 
ratio 235 U : 238 U : 239 Pu : 241 Pu :: 0.568 : 0.078 : 0.297 : 
0.057. The reactor neutrino flux is shown in Fig 2. We 
see that there although there is a significant difference 
between the reactor flux for the LMA-I and LMA-II so- 
lutions ^3 i n the geo-neutrinos energy regime. 

The Kamland experiment reports the events in (0.9 — 
1.75)MeV and (1.75 - 2.6)MeU energy. Using the the- 
oretical spectrum of the Th and U geo-neutrinos we 
have distributed the 67% of the events in the (0.9 — 
1.75)MeU bin to (0.9 - l.5)MeV and 33% of this bin 
have been assigned to the (1.5 — 2.5)MeV events. In 
the (0.9 — 1.5)MeU energy range LMA-II reactor neu- 
trinos are more suppressed than LMA-I, while in the 
(1.5 — 2.5)MeV energy range the LMA-I reactor neu- 
trinos are more suppressed than LMA-II. After this re- 
actor signal is subtracted from the KamLAND observed 
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FIG. 2: Reactor neutrino flux. Dash-dot no oscillation, 
dotted- LMA-I, continuous LMA-II. 
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FIG. 3: Simulated geo-neutrino events shown connected with 
dotted lines. Stars (LMA-I) and squares (LMA-II) denote 
observed events at KamLAND minus the reactor background. 



events, the ratio of geo events in the two energy win- 
dows , Nj/Nu will go up as we go from from LMA-I to 
LMA-II. 

In Fig 3. we plot the simulated geo-neutrino events 
shown by points connected by dotted line - where we have 
used the input values of [Th/U] = 3.8, Sin 2 29 = 0.88, 
M CC (U) = 4.2 x 10 17 kg, M oc (U) = 4.8 x 10 15 kg and 



M m (U) = 1.23xl0 17 fc<?. The geo-neutrino events have no 
dependence on Am 2 . In the same figure we have show the 
experimental points (Kamland events minus the reactor 
background) where stars represent LMA-I and squares 
represent LMA-II. Both the LMA-I and LMA-II data 
points agree with the geo-signal within the la error bars 
of the KamLAND data points. The Nj/Nu ratio which 
can be read off from the graph is considerably large for 
the LMA-II points compared to the LMA-I points which 
leads to a large value of t age for LMA-II compared to 
LMA-I. The error in determination of t age from Kam- 
LAND is estimated as follows. The systematic error in 
each of Nj and N H from Table II of 1] is 6.4% . The 
combined systematic and statistical error in Nj is 25% 
and in Njj it is 26%. The total systematic and statistical 
error in determination of Ni/Njj is 36%. The theoretical 
uncertainity in the initial r-process prediction of [Th/U] 
is 5%. The assumption that the geometrical factor Gi(E) 
is energy independent introduces an error of 1%. Added 
in quadrature the total error in t age is turns out to be 
36.5%. With this error we plot the 90% and 68.3% al- 
lowed region for t age as a function of mixing angle in 
Fig 1. We see that the LMA-II solution does not over- 
lap with the solar system age of 4.5Gyrs at 90%C.L.. 
We emphasize that the LMA-II data points fit the geo 
events in KamLAND, its the extra constraint of requir- 
ing Nj/Nu ~ 1.13 (which amounts to t age — A.dGyrs) 
which the LMA-II solution does not fulfill. We must 
add the cautionary caveat that although there is a large 
magnification in the ratio Ni/Nu in going from LMA-I 
to LMA-II, this ratio is meaningful when Nj and Nu 
are non-zero. The present KamLAND observations [l| of 
geo-neutrino events are consistent with zero at 2a, so the 
conclusions derived in this paper should be treated as la 
results. 
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